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1.  INTRODUCTION 


The  goal  of  EXCEDE  III  is  to  provide  a  validated  data  base  which  can  be  used  to  assess 
chemical  kinetic  mechanisms  in  order  to  reduce  the  uncertainties  in  the  predictions  of  the 
effects  of  nuclear  events  on  atmospheric  radiance.  Analyzed  data  from  EXCEDE  III  will  help 
identify  and  predict  optical  signatures  from  the  high-altitude  detonation  of  a  nuclear  device. 

To  this  end,  the  reaction  of  electrons  with  the  natural,  quiescent  atmosphere  resulting  from 
dosing  by  energetic  electrons  is  investigated.  A  beam  of  electrons  having  a  nominal  initial 
energy  of  2.6  keV  was  injected  into  the  upper  atmosphere  between  80  and  115  km.  Optical 
instruments  viewed  the  interaction  of  these  electrons  with  the  atmosphere  from  both  onboard 
the  electron  accelerator  module  and  from  offboard  platforms.  Measurements  were  made  of 
both  the  electron  input  and  the  resultant  atmospheric  response.  The  data  contains  information 
necessary  to  characterize  the  atmosphere  before,  during,  and  after  the  injection  of  energetic 
electrons. 

The  dose  rates  achieved  by  EXCEDE  III,  on  the  order  of  1 0"  cm  '  at  110  km,  approach 
the  dose  range  of  high  altitude  nuclear  scenarios.  The  infrared  signatures  from  this  dose 
should  be  scalable  to  nuclear  code  calculations,  though  only  a  portion  of  the  phenomenology 
is  applicable.  Current  uncertainties  in  the  phenomenological  models  based  on  reaction  rate 
coefficients  and  calculated  electron  concentrations  that  affect  infrared  emissions  are  thought  to 
lie  within  a  factor  of  2  or  3  in  the  best  known  situations  and  much  worse  in  many  other  cases. 
Although  production  rates  of  known  species  can  be  quantified  in  laboratory  experiments,  the 
total  production  of  IR  active  species  is  not  currently  available.  It  is  important,  then,  to  verify 
this  degree  of  confidence  by  providing  a  controlled  electron  dose  experiment  which  will 
provide  an  atmospheric  database  to  which  these  computer  models  can  be  compared. 

Of  particular  importance  for  system  applications  is  validation  of  both  the  emission  rate  and 
the  mechanisms  of  formation  of  important  infrared  radiators  in  the  nuclear  disturbed 
atmosphere.  Radiation  in  spectral  regions  which  normally  form  atmospheric  windows  would 
mean  that  these  windows  are  not  available  during  a  nuclear  disturbance.  The  particular 
radiating  species  has  an  important  effect  on  the  design  of  a  system.  Species  of  interest  in  this 
context  are  C02  (v2  and  v3),  NO  (Av  =  1  and  Av  =  2),  CO,  03,  and  other  oxides  of  nitrogen 
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such  as  N20  and  N02.  Emissions  from  ions,  such  as  NO+,  could  mean  geomagnetically 
formed  striations  for  a  satellite  system.  Unknown  or  disputed  reaction  rates  include  the 
quenching  of  precursor  species  such  as  N(2D)  by  electrons,  fractional  yields  of  species  such  as 
N(2P)  and  N2(A),  effects  of  these  species  in  the  disturbed  atmosphere,  and  creation 
mechanisms  for  species  such  as  highly  rotationally  excited  NO,  O-..  and  N20.  Emissions  due 
to  very  high  rotation  levels  which  are  not  quenched  in  the  field  experiment  may  have  a  great 
impact  on  the  atmospheric  bandpasses  where  current  surveillance  instruments  are  being 
designed. 

This  report  describes  the  investigation  of  NO  emission  observed  during  EXCEDE  III 
performed  under  Contract  No.  F19628-93-C-0052.  To  understand  the  NO  formation 
mechanisms,  it  is  important  to  characterize  the  formation  of  the  N  atom  precursors  thought  to 
be  the  dominant  source  of  NO.  In  particular,  the  reactions  of  N(4S),  N(2D),  and  N(2P)  with 
02  are  considered  in  some  detail.  To  this  end.  Section  2  describes  the  chemical  kinetics 
model  and  a  comparison  of  the  predicted  N(2D)  and  N(2P)  column  densities  to  those  derived 
from  the  5200  A  photometer  and  the  visible  spectrometer.  Section  3  reports  the  use  of 
quasiclassical  trajectories  to  determine  the  energy  dependent  reaction  rate  constant  and 
NO(v,J)  product  distributions  for  the  N(4S)+02  reaction,  which  has  been  proposed  as  the 
mechanism  for  rotationally  hot  NO  formation.  The  predicted  NO  populations  using  the  N 
atom  number  densities  from  the  chemical  kinetics  model  and  the  known  reaction  rate 
constants  for  the  N(2D)+02  and  N(2P)+02  reactions  and  calculated  rate  constants  described  in 
Section  2  are  compared  with  the  rotationally  thermal  and  hot  components  of  the  measured  NO 
column  densities  in  Section  3.  A  brief  summary  is  given  in  Section  4. 
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2.  ANALYSIS  OF  N(2D)  AND  N(2P)  EMISSIONS 


The  primary  source  of  NO  emissions  in  an  aurorally  enhanced  upper  atmosphere  is 
through  odd-nitrogen  chemiluminescent  reactions  involving  N(4S),  N(2D),  and  N(2P).  Recent 
papers  speculating  on  the  creation  of  the  NO  hot  bands  have  identified  high  energy  nitrogen 
atoms  as  the  source  of  these  emissions.  One  of  the  goals  of  the  EXCEDE  III  measurement 
program,  which  has  been  described  previously  is  to  understand  the  electron-beam  induced 
infrared  fluorescence  resulting  from  chemical  reactions  and  energy  transfer  processes.  To 
characterize  the  relative  importance  of  the  N  atom  precursors  in  NO  formation,  the  N(  D)  and 
N(2P)  emissions  have  been  measured  in  EXCEDE  III.  These  measurements  are  compared  to 
a  chemical  kinetics  model  which  is  tested  for  consistency  by  using  supporting  measurements, 
such  as  electron  deposition,  NO  chemiluminescence,  and  N2  second  positive  (2P)  emissions. 
The  atomic  nitrogen  formed  when  an  energetic  electron  interacts  with  the  atmosphere 
determines  not  only  how  much  NO  is  formed,  but  also  the  spectral  distribution  of  the 
chemiluminescence.  The  relative  production  of  these  atoms  in  a  electron  excited  atmosphere 
has  been  measured  under  controlled  conditions  in  EXCEDE  III.  Simultaneous  measurements 
of  the  2P  — ►  4S,  2D  — ►  4S,  and  N2  first  negative  (IN)  bands  have  been  performed  during  this 
experiment.  A  carefully  controlled  electron  beam  dosing  profile  gave  a  direct  measure  of  the 
odd  nitrogen  production  rates  as  a  function  of  electron  injection.  The  energy  deposition 
profile  was  measured  using  3914  A  emission,  and  is  well  represented  by  a  generalized 
Gaussian  function/2^  Since  the  exact  mission  profile  was  known,  the  emissions  within  the 
field-of-view  of  each  detector  could  be  compared  with  the  chemical  kinetics  responsible  for 
the  nitrogen  formation.  This  report  presents  an  analysis  of  the  nitrogen  atoms  produced  in  its 
two  lowest  metastable  excited  states  and  predicts  the  amount  of  N(4S)  formed  based  on 
branching  ratios  reported  in  the  literature. 

The  instruments  used  to  obtain  the  N  atom  concentrations  consisted  of  a  visible 
spectrometer  for  N(2P),  and  dedicated  photometer  for  N(2D).  A  second  photometer,  looking 
at  N2  IN,  was  boresighted  with  each  of  these  instruments.  Since  several  different  instruments 
were  used  for  these  measurements,  they  are  described  separately  in  the  sections  below.  In  the 
following  sections,  we  present  the  measurement  of  N(2D)  and  N(2P),  the  modeled 
concentrations  based  on  the  observed  N2  (0,0)  IN  emissions,  and  finally  a  discussion  on  the 
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comparison. 


2.1  Measurement  of  Metastable  N(2D)  and  N(2P) 

2.1.1  N(2D) 

The  N(2D)  emission  was  measured  by  a  filter  photometer  with  a  center  wavelength  at 
5200  A  and  a  full  width  at  half  maximum  (FWHM)  of  6.06  A.  The  field-of-view  (FOV)  of 
the  photometer  was  nominally  6°  and  was  pointed  18°  into  the  afterglow.  Data  were  collected 
at  a  sample  rate  of  750  Hz.  In  order  to  correct  for  possible  interfering  emission,  such  as  the 
0-3  band  of  the  N2  IN  series  at  5228  A  or  the  Herzberg  bands  of  02,  an  additional 
photometer  with  a  center  wavelength  at  5228  A  with  a  FWHM  of  5.79  A  was  co-aligned  with 
the  5200  A  photometer. 

2.1.2  N(2P) 

The  measurement  of  N(2P)  produced  in  EXCEDE  III  was  made  by  measuring  its  relative 
intensity  to  known  features  using  a  visible  spectrometer.  A  more  detailed  description  of  this 
instrument  is  reported  by  Reider  et  al .  ■ 2 1  This  instrument  had  a  spectral  scan  time  of  2.8  s. 
Since  alternating  scans  used  a  neutral  density  filter  (OD  =  2)  to  extend  the  instrument 
dynamic  range,  one  3466  A  measurement  was  made  every  5.6  s.  This  time,  convolved  with 
beam-on  times  of  7.1  s,  limited  the  number  of  N(2P)  scans  to  15,  6  on  upleg,  one  at  apogee, 
and  8  on  the  downleg.  The  feature  at  3466  A  is  in  the  region  of  the  spectrum  where  the  N2 
and  N2  2nd  Positive  bands  are  located.  In  fact,  the  Av  =  1  series  of  N2(2+)  has  a  band  in 
near  coincidence  with  N(2P),  at  3468  A. 

2.2  Chemical  Kinetics  Model 

The  most  important  infrared  emitters  in  an  electron-dosed  atmosphere  include  NO,  C02, 
NO+,  CO,  and  03.  A  general  chemical  kinetics  model  has  been  developed  to  describe  the 
chemical  and  collisional  processes  occurring  during  the  dosing  of  the  atmosphere  with  an 
energetic  electron  beam.  The  model  uses  a  finite  rate  approach  to  the  modeling  of  the 
generation  of  energy  dependent  secondary  electrons,  vibrationally  and  electronically  excited 
species  resulting  from  collisional  energy  transfer  and  chemiluminescent  reactions.  In 
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particular,  time  and  energy  dependent  rate  equations  are  solved  for  arbitrary  electron  beam 
dosing  to  obtain  species  number  densities  as  a  function  of  time.  The  secondary  electron 
spectrum  is  treated  by  setting  up  an  energy  grid  (histogram)  at  sufficiently  high  resolution  to 
accurately  solve  electron-molecule  and  electron-electron  interactions.  Time-dependent  rate 
equations  are  then  solved  for  each  energy  bin.  Chemical  and  collisional  processes  for 
individual  vibrational  states  are  treated  explicitly.  Self-relaxation  of  secondary  electrons, 
which  is  important  due  to  the  high  electron  beam  fluxes  produced  in  the  EXCEDE 
experiment,  has  been  treated  by  solving  the  time-dependent  Fokker-Planck  equation.*3,4*  In 
the  remainder  of  this  section,  the  detailed  chemical  and  collisional  processes  influencing  the 
production  and  destruction  of  metastable  nitrogen  atoms  is  discussed  in  some  detail. 

When  electrons  interact  with  the  natural  atmosphere,  dissociation  and  ionization  form  odd 
nitrogen  in  the  form  of  N(2P),  N(2D),  and  N(4S).  Collisions  between  these  species  and  02 
create  NO.*5*  This  production  rate  is  in  turn  found  to  be  dependent  on  the  ratio  of  N(2D)  to 
N(4S)  formed  in  the  interaction  of  electrons  and  N2.  Nitric  oxide  is  currently  thought  to  be 
produced  primarily  by  the  reaction  of  N(2D)  with  02 

N(2D)  +  02->  NO  +  O  (1) 

where  the  reaction  rate  constant*6*  and  NO  vibrational  distribution*7,8*  are  well  established  at 
room  temperature.  Recent  laboratory  studies*8*  have  indicated  that  N(2P)  produces  NO  in 
high  vibrational  and  rotational  states 

N(2P)  +  02-+  NO  +  O  (2) 


with  a  rate  constant  a  factor  of  2  smaller  than  Reaction  (l).*6*  Furthermore,  it  has  been 
suggested  by  Sharma  et  al.,*9)  that  translationally  hot  N(4S)  atoms*10*  may  be  a  significant 
source  of  rotationally  hot  NO  in  the  thermosphere 
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(3) 


N(4S )*  +  02~+  NO  +  O 


Although  the  reaction  rate  constant  and  the  final  vibrational-rotational  distribution  at  high 
translational  energies  have  not  been  measured,  recent  calculations  by  Duff  et  al.,(11)  have 
established  realistic  rate  constants  and  NO  vibrational-rotational  distributions  for  Reaction  (3). 
A  realistic  description  of  the  formation  and  emission  of  NO  in  the  ambient  or  auroral 
atmosphere  requires  a  detailed  treatment  of  nitrogen  atom  chemistry. 

Nitrogen  atoms  are  produced  by  several  major  processes;  dissociation  of  N2  by  the 
primary  electrons  (ep)  in  the  electron  beam,  Reaction  (4),  dissociation  of  N2  by  the  secondary 
electrons  (es)  produced  by  ionization  of  atmospheric  species  by  the  electron  beam,  Reaction 
(5),  and  recombination  reactions  of  the  primary  ionic  products  with  secondary  electrons. 
Reactions  (6)  and  (7), 


e  +  N  -+  e  +  N(S)  +  N(2D)  (4a) 

pip 

e  +  N  -*  e  +  N(4S)  +  N(2P)  (4b) 

pip 

e  +  N  ->  e  +  e  +  N+  +  ViN(S)  +  lAN(2D)  (4c) 

pips 

e  +  N  -*  e  '  +  N(4S)  +  N(2D )  (5a) 

s  2  s 

e  +  N  -*■  e  '  +  N(4S)  +  N(2P)  (5b) 

s  2  s 
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(6) 


N~  +  e 

L  s 


(2  -  ax  )N(4S)  +  aN(2D) 


NO*  +  e  -»  (1  -  a  )N(4S)  +  a  N(2D )  +  O 

s  2  2 


(7) 


The  total  production  rate  of  N  atoms  from  Reactions  (4a)  and  (4b)  is  assumed  to  be 
1.5P(N2),(I2’  ,3)  where  P(N^)  is  the  production  rate,  with  the  relative  formation  rate  of 
N(2D)  and  N(2P)  in  the  ratio  0.52:0.48.(14^  For  dissociative  ionization  by  primary  electrons. 
Reaction  (4c),  the  production  rate  is  0.50P(N2).(15)  Thus,  the  total  production  rate  of  N  atoms 
from  Reaction  (4)  is  2PCN2),  or  1 .24  N  atoms/ion  pair,  in  good  agreement  with  the  results  of 
Porter  et  al./12)  and  Rusch  et  al.,(13)  The  cross  sections  for  dissociation  of  N2  by  low  energy 
secondary  electrons.  Reaction  (5),  are  based  on  the  analysis  of  Zipf  et  al./14^  The  rate 
constants  for  the  recombination  of  N2  and  NO+  ions  with  secondary  electrons,  Reactions  (6) 
and  (7),  have  been  measured  by  Mehr  and  Biondi*16^  and  Alge  et  al./17^  respectively.  The 
associated  branching  ratios  for  the  Nj  and  NO+  recombination  reactions  are  a,  =  1.86(18)  and 
a2  =  0.76.<l9)  In  addition  to  Reactions  (4)  through  (7),  minor  sources  of  N(2D)  considered  are 
the  charge  exchange  reactions  of  N2  with  0(2°^  and  N+  with  02^21'22^ 


N2  +° 


N(D)  +  NO 


(8) 


N*  +  02->  N(  D)  +  02 


(9) 


Nitrogen  atoms  are  also  subject  to  quenching  by  other  species  which  compete  with 
Reactions  (1)  and  (2)  in  the  production  of  NO 


N(2D)  +  N  -*>  N(4S)  +  N2 


(10) 
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(11) 


N(2D )  +  O  —■  N(4S)  +  o 

N(2p )  +  o -*■  N(2D )  +  O  (l2a) 

N(2p )  +  O  NO*  +  e  .  (12b) 

The  quenching  rates  of  Reactions  (10)  -  (12)(6,23)  are  slow  on  the  time  scale  of  the  EXCEDE 
III  measurements  in  comparison  to  the  efficient  quenching  of  N(  D)  and  N(  P)  by  secondary 
electrons, (24)  which  results  from  the  large  density  of  secondary  electrons  created  by  the 
ionization  rates  achieved  in  the  experiment 

N(2D)  +  e  -*•  N(*S)+e  (13) 

S 

N(2P)  +  e  ->  N(4S,2D )  +  e  '  (14) 

s  s 


In  addition  to  collisional  quenching,  N(2D)  and  N(2P)  can  undergo  radiative  relaxation 

N(2D )  -*■  N(4S)  +  h\  (I5) 

N(2P)  ->  N(2D,4S)  +  hv  .  (I6) 

Since  the  N  atom  production  from  direct  electron  impact  excitation  and  recombination 
reactions  depends  on  the  low  energy  secondary  electron  spectra,  it  is  important  that  the 
histogram  treatment  retain  the  major  features  of  the  secondary  electron  spectra.  Secondary 
electrons  are  formed  via  collisional  ionization  of  M(=  N2,  02,  or  O) 
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(17) 


e  +  M  -*  e  +  e  +  products 

p  p  * 

The  nascent  secondary  electron  energy  distribution  resulting  from  (17)  is  that  given  by  Opal 
et  al./25*  which  then  relaxes  by  collisions  with  major  species 

e  +  M  —*  e  1  +  M  (18) 

S  S 


and 


e  +  M  -*  e  '  +  e  "  +  M+  (19) 

s  s  s 

as  well  as  electron-electron  collisions 

e  +  e  '  -*■  e  "  +  e  (20) 

s  s  s  s  v  7 

where  the  primes  refer  to  secondary  electrons  in  different  energy  bins.  The  ion  production 
rates  from  Reaction  (17)  are  related  to  the  ion  pair  production  rate  via  the  expressions  given 
by  Rees  and  Jones. ^26)  The  rate  constants  for  Reactions  (18)  and  (19)  were  taken  from  the 
reviews  of  Itikawa  et  al.,(27,28)  for  N2  and  02  and  Itikawa  and  Ichimura<29)  for  O. 

One  of  the  most  important  inputs  to  the  model  is  the  electron  beam  characteristics,  i.e.,  the 
radial  profile  of  the  ion-pair  production  rate.  Historically,  the  electron  deposition  information 
has  been  deduced  from  the  N2  first  negative  3914  A  emission,  where  it  is  assumed  that  there 
are  14.1  N2  ions  created  for  each  3914  A  photon/26^  Therefore,  the  beam  ion  pair 
production  rate  profiles  used  in  the  current  model  are  based  on  an  analysis  of  the  3914  A  data 
by  Rieder  et  al.,(2)  The  N2  production  rate,  P(N2),  is  related  to  the  total  ion  pair  production 
rate,  Pion,  by 


0.92  [Ay 

1.15[Ay  +  1.5 [<9]  +  0.56[O]  P‘on 


(21) 
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where  the  numerical  factors  in  the  expression  represent  the  relative  ionization  cross  sections 
for  N2,  02,  and  0.(26)  Using  the  atmospheric  profile  appropriate  for  the  EXCEDE  III 
experiment,  the  relative  production  rate  of  N2  ions  is  approximately  independent  of  altitude 
(ranging  from  0.60  at  95  km  to  0.62  at  1 15  km). 

Representative  predictions  of  the  production  and  destruction  rates  for  N(  D)  and  N(  P)  are 
presented  in  Figures  1  through  4  for  an  altitude  of  1 15  km.  The  electron  beam  ion  pair 
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production  rate  is  represented  by  a  generalized  Gaussian  function(2)  with  a  peak  ion  pair 
production  rate  of  1.6  x  10n  ion  pairs/cm3-s  and  a  full  width  at  half  maximum  of  14  ms. 


-0.1  0.0  0.1  0.2  0.3  0.4  0.5 


TIME  (S) 

Figure  1.  Major  Processes  for  the  Production  Rate  of  N(2D)  at  115  km. 
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TIME  (S) 

Figure  2.  Major  Processes  for  the  Destruction  Rate  of  N(2D)  at  115  km. 
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2.3  Comparison  of  Data  with  The  Kinetics  Model 


The  chemical  kinetics  mechanism  described  previously  has  been  incorporated  into  the 
EXCEDE  chemical  kinetics  data  base.  The  resultant  data  base  contains  6703  reactions  and 
energy  transfer  processes.  The  chemical  kinetics  equations  are  integrated  using  Gear’s  method 
for  stiff  differential  equations^30*  as  a  function  of  time  for  a  specified  electron  beam  ion-pair 
production  rate  and  altitude.  The  atmospheric  profiles  for  N2,  02,  and  O  were  obtained  from 
the  MSISE-90  atmospheric  model(31)  for  the  geophysical  parameters  appropriate  to  the 
EXCEDE  III  experiment. 

The  geometry  of  the  flight  experiment,  which  includes  the  trajectory  of  the  gun  module, 
the  geomagnetic  azimuth  and  dip  angles,  and  the  viewing  geometry  of  the  3466  A  visible 
spectrometer  and  5200  A  photometer  on  the  sensor  module,  has  also  been  incorporated  in  the 
model.  The  time-dependent  nitrogen  atom  number  densities  are  integrated  along  the 
instrument  line-of-sight  to  obtain  the  column  densities  as  a  function  of  time.  The  resultant 
column  densities  are  then  integrated  over  the  6°  FOV  of  the  instruments  for  comparisons  with 
the  column  densities  obtained  from  3466  A  and  5200  A  data. 

To  illustrate  the  electron  beam  energy  deposition  throughout  the  EXCEDE  flight,  the  total 
number  of  ion  pairs  produced  by  the  electron  beam  is  shown  as  a  function  of  mission  elapsed 
time  (or  altitude)  in  Figure  5.  The  number  of  ion  pairs  was  obtained  by  integrating  the 
observed  3914  A  scanning  photometer  data  along  the  component  of  gun  vehicle  velocity 
vector  perpendicular  to  the  magnetic  field  lines.  The  maximum  number  of  ion  pairs  are 
created  near  102  km  on  upleg,  where,  due  to  the  near  zero  velocity  across  the  field  lines,  the 
atmosphere  is  dosed  for  nearly  0.5  seconds.  The  average  electron-deposition  time,  which  is 
simply  the  root  mean  square  electron  beam  diameter  divided  by  the  velocity  across  the  field 
lines,  is  shown  in  Figure  6.  As  altitude  increases,  the  number  of  ion  pairs  produced  decreases 
due  to  the  shorter  deposition  times  (i.e.,  the  velocity  across  the  field  lines  is  increasing  with 
altitude),  see  Figure  6.  On  downleg,  the  number  of  ion  pairs  again  increases  as  altitude 
decreases,  although  the  total  number  of  ion  pairs  created  is  substantially  less  on  downleg  than 
on  upleg.  This  effect  is  due  to  the  atmosphere  being  irradiated  by  the  electron  beam  for  times 
up  to  a  factor  of  50  longer  on  upleg  than  on  downleg,  again  as  shown  in  Figure  6.  The  two 
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data  points  at  approximately  114  and  113  km  are  about  75%  of  the  expected  beam  power 
output  due  to  beam  load  faults. 

The  comparisons  of  the  N(2D)  and  N(2P)  column  densities  obtained  from  the  5200  A  and 
3466  A  data  with  the  model  predictions  which  have  been  integrated  along  the  line-of-sight 
and  integrated  over  the  sensor  field-of-view  are  shown  in  Figures  7  and  8,  respectively.  The 
agreement  between  the  model  predictions  and  data  is  excellent  on  both  upleg  and  downleg  for 
the  N(2D)  measurements.  The  altitude  dependence  of  the  N(2D)  column  density  measurement, 
which  looks  1 8°  into  the  afterglow,  is  controlled  by  quenching  of  N(2D)  by  secondary 
electrons  with  increasing  importance  of  the  N(2D)+02  reaction  near  100  km.  On  upleg,  the 
largest  discrepancy  between  the  model  and  data  is  at  109  and  111  km  on  upleg,  which  is  not 
understood,  although  it  is  likely  related  to  the  uncertainty  in  the  secondary  electron  quenching 
rate  of  N(2D).  The  differences  between  upleg  and  downleg  N(2D)  column  densities  simply 
reflect  the  differences  the  ion  pair  production  rate,  although  the  ratio  of  upleg  to  downleg 
N(2D)  column  densities  is  less  sensitive  to  ion  pair  production  rate  due  to  the  significant 
quenching  of  N(2D)  and  the  fact  that  the  primary  production  process  for  N(2D)  is 
recombination  of  secondary  electrons  with  N2.  The  production  rate  of  N(2D)  from 
recombination  roughly  increases  as  the  square  root  of  the  ion  pair  production  rate. 

The  comparison  of  the  model  predictions  with  the  3466  A  data  shown  in  Figure  8  indicate 
that  the  N(2P)  chemistry  is  not  as  well  established  as  the  N(2D)  chemistry.  This  comparison 
would  indicate  that  the  assumed  production  rate  of  0.22  N(2P)/ion  pair  used  in  Reaction  (4b) 
is  approximately  a  factor  of  2  low.  Since  the  total  production  rate  of  N  atoms  is  fairly  well 
established,  increasing  the  rate  for  Reaction  (4b)  by  a  factor  of  2  would  imply  that  only  N(  P) 
is  produced  by  dissociation  of  N2  by  primary  electrons.  Such  a  modification  in  the 
N(2D)/N(2P)  relative  production  rate  from  Reaction  (4)  would  not  significantly  alter  the 
excellent  agreement  between  the  predictions  and  the  5200  A  data  since  the  major  production 
mechanism  for  N(2D)  is  the  recombination  of  N2  with  secondary  electrons,  Reaction  (6),  as 
shown  in  Figure  1.  Furthermore,  the  increase  in  N(2P)  would  provide  more  N(2D)  due  to  the 
increased  quenching  rate  from  Reaction  (14). 
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Figure  5. 


Figure  6. 


Total  Number  of  Ion  Pairs  Along  the  Line-of-Sight  as  a  Function  of  Mission 
Elapsed  Time  for  the  EXCEDE  III  Experiment. 


DOSE  TIME  (S) 

Calculated  Dose  Time  as  a  Function  of  Mission  Elapsed  Time  for  the  EXCEDE 
III  Experiment. 
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Figure  7. 


Figure  8. 
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Comparison  of  the  Time  Dependence  of  the  Calculated  N(2D)  Column  Density  (o) 
with  the  5200  A  Photometer  Data  (•). 


Comparison  of  the  Calculated  Time  Dependence  of  the  N(2P)  Column  Density  (o) 
with  the  3466  A  Spectrometer  Data  (•). 
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2.4  Summary 


The  EXCEDE  III  experiment  has  provided  quantitative  characterization  of  the  N(2D)  5200 
A  and  N(2P)  3466  A  emissions  resulting  from  the  irradiation  of  the  lower  thermosphere  (90  - 
115  km)  by  an  intense  electron  beam.  A  detailed  chemical  kinetics  model  indicates  that  the 
N(2D)  emission  can  be  quantitatively  explained  by  assuming  accepted  production  rates  of 
N(2D)  and  N2  created  by  the  electron  beam  dissociation,  consistent  with  previous  auroral 
models.  Using  the  previous  estimate  of  0.22  N(2P)  per  ion  pair*14*  from  electron  dissociation 
of  N2  results  in  a  model  prediction  which  is  approximately  50%  below  the  data,  in 
disagreement  with  previous  analyses  of  N(2P)  auroral  and  dayglow  emission.*14’32,33* 

However,  it  should  be  recalled  that  the  EXCEDE  III  measurement  of  N(2P)  was  not  made 
under  steady  state  conditions,  which  are  appropriate  to  the  atmospheric  observations.  The 
present  analysis  would  imply  that  the  amount  of  N(2P)  from  N2  dissociation  should  be  about  a 
factor  of  2  greater  than  previously  assumed  or  an  additional  source  of  N(2P)  exists  for  intense 
electron  beams  such  as  that  used  in  EXCEDE  III. 
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3.  QUASICLASSICAL  TRAJECTORY  STUDY  OF  THE 
TRANSLATIONAL  ENERGY  DEPENDENCE  OF  THE 
N(4S)  +  02(X3S*)  -*•  NO(X2II)  +  0(3P)  REACTION 


3.1  Introduction 

It  is  recognized  that  nitric  oxide  plays  an  important  role  in  the  energy  budget  and 
photochemistry  of  the  lower  thermosphere.  The  reactions  of  metastable  nitrogen  atoms  with 
molecular  oxygen  are  thought  to  be  the  major  source  of  NO  chemiluminescence  in  the 
thermosphere/5*  Standard  kinetic  models  of  NO  formation  assumed  that  the  thermal 
chemistry  of  N(2D)  is  the  most  significant  contributor  to  NO  formation,  while  the  relaxation 
of  N(2P)  to  N(2D)  and  the  reaction  of  N(2P)  with  02  make  a  relatively  minor  contribution  to 
NO  formation  in  the  thermospheric  dayglow(33)  and  aurora/14’32*  At  thermospheric 
temperatures,  the  reaction  of  N(4S)  with  02  is  much  slower  than  Reactions  (1)  and  (2)  and, 
therefore,  has  not  been  considered  a  significant  source  of  NO.  However,  Solomon*34*  has 
suggested  that  the  reaction  of  translationally  hot  N(4S)  atoms  with  02  may  be  an  important 
factor  in  the  formation  of  NO  in  the  thermosphere.  Detailed  models  involving  the  solution  of 
the  Boltzmann  equation  for  the  kinetic  description  of  N(4S)  atoms  in  the  thermosphere  have 
shown  the  translational  energy  distribution  to  be  highly  nonthermal/10’  35'37) 

Model  predictions  by  Gerard  and  coworkers*10’36*  have  also  shown  that  the  reaction  of 
these  nonthermal  N(4S)  with  02  makes  a  significant  contribution  to  thermospheric  NO 
formation  in  the  daytime.  Furthermore,  Sharma  et  al.,(9)  have  argued  for  the  importance  of 
the  N(4S)  +  02  reaction  in  recent  observations  of  highly  rotationally  excited  NO  vibrational 
emissions  in  the  dayglow.*38’ 39)  Their  analysis  was  based  on  two  recent  studies;  a 
quasiclassical  trajectory  study  of  ground  state  N(4S)  +  02  reaction  by  Gilibert  et  al.,*40* 
indicating  that  the  N(4S)  reaction  efficiently  produces  NO  in  high  ro-vibrational  states  and  the 
calculated  nonthermal  thermospheric  N(4S)  translational  energy  distribution/ 10)  Further 
evidence  showing  the  importance  of  Reaction  (3)  in  producing  highly  rotationally  excited  NO 
in  the  atmosphere  has  been  provided  by  the  analysis  of  NO  emission  observed  during  an 
artificial  auroral  experiment  in  the  lower  thermosphere/41*  However,  Rawlins  et  al.,(8)  have 
interpreted  the  formation  of  highly  rotationally  excited  NO  in  laboratory  experiments  from  the 
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reaction  N(2D),  although  it  has  been  pointed  out  by  Smith  and  Ahmadjian(38)  that  the  density 
of  N(2P)  is  insufficient  to  account  for  the  atmospheric  observations.  In  order  to  assess  the 
importance  of  the  N(4S)  +  02  as  a  source  of  vibrationally  and  rotationally  excited  NO,  the 
reaction  rate  constant  and  NO  vibrational  -rotational  distribution  must  be  determined  as  a 
function  of  N(4S)  translational  energy. 

The  present  theoretical  study  uses  extensive  quasiclassical  trajectory(QCT)  calculations  to 
determine  the  reaction  rate  constant  and  the  nascent  NO  vibrational-rotational  distributions 
formed  by  the  translationally  hot  N(4S)  +  02  reaction  as  a  function  of  N(4S)  translational 
energy  from  0.5  to  2.75  eV.  A  new  analytical  fit  to  the  2A'  and  4A'  ab  initio  potential  energy 
surfaces  of  Walch  and  Jaffe(42^  is  used  to  predict  the  reaction  attributes  of  N(4S)  +  02.  The 
reaction  attributes  have  been  independently  calculated  on  each  potential  energy  surface  and 
then  statistically  averaged  using  the  electronic  degeneracy  factors.  The  lowest  A'  potential 
energy  surface,  which  has  an  estimated  activation  energy  of  approximately  0.3  eV,(40)  makes 
the  dominant  contribution  to  the  reaction  rate  for  temperatures  below  1500  K.  Although  the 
4 A'  surface,  which  has  an  estimated  barrier  of  0.65  eV,(42)  does  not  contribute  to  the  reaction 
rate  for  temperatures  below  1 500  K,  it  must  be  considered  for  the  translationally  "hot"  N(4S) 
reaction.  The  results  of  this  study  provide  the  reaction  attributes  of  N(4S)  +  02  useful  for 
modeling  the  production  of  NO  from  "hot"  N(4S)(10’  35‘37)  and  the  NO  day  glow  emission 
spectra/9,  38,  39) 

The  only  experimental  data  on  the  N(4S)  +  02  reaction  are  thermal  rate  constant 
measurements, (43)  which  are  reliable  up  to  1500  K,  the  NO  vibrational  distribution  at  room 
temperature/44'46^  and  a  preliminary  study  of  energy  dependence  of  the  NO  vibrational- 
rotational  emission  using  bandpass  filters  covering  the  NO  fundamental  and  overtone  spectral 
regions/47^  No  direct  experimental  measurements  of  the  energy  dependence  of  the  reaction 
cross  section,  or  equally  important  for  thermospheric  modeling  of  NO  emission  spectra,  the 
disposal  of  available  energy  into  NO  vibration  and/or  rotation,  are  currently  available.  There 
also  is  evidence  that  the  measured  NO  vibrational  distributions  have  undergone  vibrational 
relaxation,  and  therefore  may  be  colder  than  the  nascent  distribution. 
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There  have  been  two  previous  quasiclassical  trajectory  calculations  reported  for  the  N(4S) 

+  02  reaction.  The  first  study  was  performed  by  Jaffe  et  al.,(48)  using  the  2  A'  and  4A' 
surfaces  mentioned  above.  However,  their  study  emphasizes  the  thermal  reaction  rate  constant 
at  high  temperatures  and,  although  consistent  with  experiment,  do  not  provide  rate  constants 
as  a  function  of  collision  energy  or  final  state  distributions  of  the  products.  The  second 
quasiclassical  trajectory  study  by  Gilibert  et  al.,(l3)  only  considered  the  lowest  potential  energy 
surface,  2A',  and  therefore  is  applicable  to  thermal  collisions  below  1500  K  or  collision 
energies  less  than  approximately  0.75  eV.  More  importantly,  the  fitting  of  an  analytical 
functional  form  to  the  calculated  ab  initio  2 A'  potential  energy  surface  introduced  an  artificial 
barrier  of  approximately  0.55  eV  in  the  reactant  channel,  which  implies  that  the  calculated 
threshold  for  reaction  is  too  high  by  approximately  0.22  eV.  Thermal  rate  constants  calculated 
using  the  2 A'  PES  of  Gilibert  et  al.,  would  result  in  an  activation  energy  approximately  5 
kcal/mole  higher  than  experiment.  Recently,  Gilibert  et  al./49*  have  compared  their  QCT 
results  with  a  Reactive  Infinite  Order  Sudden  Approximation  (R-IOSA)  study  of  the  N(4S)  + 
02  reaction.  The  general  features  of  the  reaction  dynamics  obtained  from  the  R-IOSA  study 
are  in  general  agreement  with  the  QCT  calculations. 

The  remainder  of  this  section  is  organized  as  follows.  Section  3.2  discusses  the 
construction  of  the  potential  energy  surfaces  used  to  study  the  reaction  of  hyperthermal  N(4S) 
and  with  02  in  producing  NO.  The  trajectory  calculations  are  briefly  discussed  in  Section 
3.3,  while  in  Section  3.4  results  of  the  chemical  dynamic  study  are  presented.  A  brief 
summary  is  given  in  Section  3.5. 

3.2  Potential  Energy  Surfaces 

The  reaction  of  N(4S)  +  02(X  3S‘)  to  form  NO(X2II)  +  0(3P)  is  exothermic  by  1.38  eV 
and  occurs  on  two  electronic  potential  energy  surfaces  of  2 A'  and  4A'  symmetry  (neglecting 
spin-orbit  coupling).  The  electronic  degeneracies  are  such  that  4/12  of  the  N(4S)  +  02 
collisions  occur  on  the  4A'  surface  compared  to  2/12  on  the  ground  state  2 A'  surface.  The 
remaining  6/12  of  the  collisions  occur  on  the  6A'  surface  which  correlates  with  the  excited 
state  products  NO(a4II)  +  0(3P),  and  is  not  considered  further.  An  extensive  set  of  ab  initio 
quantum  mechanical  calculations  for  the  2 A'  and  4A'  potential  energy  surfaces  have  been 


20 


carried  out  by  Walch  and  Jaffe*42*  to  characterize  the  saddle  point  geometries  and  to  a  lesser 
extent  the  minimum  energy  path.  Other  previous  work  on  the  N02  system  has  been  discussed 
in  some  detail  by  Gilibert  et  al.,*40* 

The  recent  ab  initio  calculations  of  Walch  and  Jaffe  for  the  2 A'  and  4A'  potential  energy 
surfaces  (PES)  are  used  as  a  basis  for  the  analytical  representation  of  the  N(4S)  +  02  reaction. 
Their  study  involved  a  total  of  66  ab  initio  points  for  the  A'  surface  and  23  for  the  A' 
surface  at  the  complete  active  space  self-consistent  field  (CASSCF)  level  followed  by 
multireference  contracted  configuration  interaction  (CCI)  calculations.  An  additional  1 1 
points  were  computed  at  the  ext.CCI  level  to  approximately  characterize  the  minimum  energy 
path  and  energy  release  on  the  2 A'  ground  state  surface.  The  "best"  theoretical  estimate  of  the 
barrier  height  for  the  2 A'  surface  was  felt  to  be  2  -  3  kcal/mole  too  high  compared  with  the 
experimental  activation  energy.  In  creating  the  analytical  representation  of  the  PES,  the 
calculated  ab  initio  points  in  the  vicinity  of  the  saddle  point  were  adjusted  by  a  constant 
factor  such  that  the  barrier  heights  for  the  2 A'  and  4A'  are  in  agreement  with  the  experimental 
estimate  for  the  2A'  surface*43*  and  theoretical  estimate  for  the  4A'  surface.*42*  The  adjusted 
ab  initio  points  were  then  fit  to  an  analytical  representation  of  the  potential  energy  surface 
suggested  by  Sorbie  and  Murrell(SM)*50)  using  a  nonlinear  least  squares  method. 

The  SM  analytical  form  for  the  PES  is  given  by  the  sum  of  two-body  diatomic  potentials, 
V*2)(r),  and  a  three-body  term,  Vj3*’ 


'W'V'vV 


(2)  (2)  (2)  (3) 

W  +  FO0(r 2)  +  W  +  V  1  (ri,r2’r3) 


(22) 


where  r,,  r2,  r3  are  the  NO,  02,  NO  intemuclear  distances,  respectively.  The  diatomic 
potentials  are  represented  by  a  modified  Rydberg  potential*51* 
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where  De  is  the  dissociation  energy,  re  is  the  equilibrium  intemuclear  distance,  and  (a,,  a2,  a3) 
are  parameters  fit  to  the  RKR  potential.  The  parameters  for  NO  and  02  given  in  Table  1  are 
identical  to  those  used  by  Gilibert  et  al.,(40)  The  three-body  term  is  a  product  of  a  range 
function, 
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Table  1.  Parameters  for  The  Diatomic  extended-Rydberg  Potentials.2 


De 

re 

al 

a2 

a3 

NO(X2III/2) 

6.6144 

1.1508 

5.035 

5.151 

2.998 

5.2132 

1.2075 

5.476 

7.950 

6.341 

aEnergies  are  given  in  eV  and  distances  in  Angstroms. 


which  ensure  the  correct  asymptotic  behavior,  and  a  polynomial  expressed  as  a  function  of 
internal  coordinates 


V°\r,r,r)  =  V  T  T  T  (l  +  Tc  s'sV ) 

v  l’  2’  y  0  1  2  3V  u  jkl  1  2  3  7 


k  /, 


(24) 


where  the  internal  coordinates  are  defined 

5  =  l/v/2(pi  +  p3)  (25a) 

*2  =  P2  (25b) 

=  l/v/2(pi  -  p3)  (25c) 


with 


p  =  r  -  r :  .  (26) 

i  i  l 

The  sum  in  Equation  (24)  consists  of  a  polynomial  in  internal  coordinates  such  that  0  < j  +  k  + 
1  <  M,  where  M  =  4  for  the  present  case.  The  range  functions,  which  decay  exponentially  as  s, 
—*  co,  are  defined  in  terms  of  hyperbolic  tangents 
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(27) 


T  =  1  -  tanh(ys/2) 

I  i  l 


The  three-body  potential  defined  by  Equations  (24)  -  (27)  requires  the  determination  of  41 
parameters,  14  of  which  vanish  due  to  symmetry,  and  3  of  which  are  specified  (r°,  r2,  r°  =  r°). 
The  remaining  24  parameters  are  determined  by  a  nonlinear  fit  to  the  adjusted  ab  initio  points. 

Using  the  original  77  ab  initio  points  for  the  2 A'  resulted  in  an  analytical  function,  which 
displayed  spurious  wells  and  barriers  in  the  entrance  channel  similar  to  those  found  in  the  fit  by 
Gilibert  et  al.,(40)  This  undesired  result  was  remedied  by  generating  additional  points  in  the  N(4S) 
+  02  entrance  channel  assuming  that  the  three  body  interaction  term  decreases  exponentially  as 
rN0  — *•  oo5  This  assumption  has  no  influence  on  the  reaction  dynamics,  but  removes  the  spurious 
artifacts  in  the  analytical  fit.  The  parameters  (Yj  and  cjk|)  in  Equations  (24)  and  (27)  are 
determined  by  fitting  to  the  set  of  adjusted  ab  initio  points,  augmented  with  additional  points  in 
the  entrance  channel,  with  a  standard  deviation  of  0.047  eV  for  both  surfaces.  The  final 
parameters  for  the  2 A'  and  4A'  three-body  potential  energy  surfaces  are  given  in  Table  2. 

Contour  plots  of  the  2 A'  and  4A'  potential  energy  surfaces  are  shown  in  Figures  9  and  10  for 
a  bond  angle,  0NOO,  of  110°.  The  resultant  barrier  heights  of  the  analytical  functions  are  0.3  eV 
and  0.65  eV  for  the  2 A'  and  4A'  potential  energy  surfaces,  respectively.  The  saddle  point 
geometry  is  given  by  rN0  =  1.921  A,  r00  =  1.216  A,  and  0NOq  =  111°  f°r  the  2 A'  surface  and 
rN0  =  1.813  A,  Tqq  =  1.232  A,  and  0NOO  =  107°  for  the  4A',  in  excellent  agreement  with  the  ab 
initio  results.  The  angular  dependence  of  the  barrier  for  both  surfaces  is  also  accurately 
reproduced  by  the  analytical  function,  as  is  shown  in  Figure  1 1 .  Both  surfaces  are  characterized 
by  attractive  energy  release,(52)  which  is  typical  of  surfaces  with  early  barriers  in  the  entrance 
channel.  However,  it  is  apparent  from  the  contour  plots  that  the  curvature  of  the  4A'  minimum 
energy  path  is  smaller  than  that  for  the  2A'  which  would  imply  that  NO  vibrational  excitation 
would  be  more  efficient  for  the  2A'  PES.(53)  It  should  be  noted  that  although  the  ab  initio 
calculations  indicate  that  the  curvature  of  the  minimum  energy  path  for  4A'  is  smaller  than  that 
for  the  2A',  the  location  and  curvature  of  the  minimum  energy  path  and  the  energy  release  were 
not  completely  determined  for  the  4A'  in  the  study  of  Walch  and  Jaffe. 
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Table  2.  Three  Body  Parameters  for  the  N02  2A"  and  4A"  Potential  Energy  Surfaces 


1.701 

1.815 

-0.730 

-0.372 

3.632 

3.751 

0.0 

0.0 

1.415 

1.239 

-2.252 

-1.627 

0.0 

0.0 

13.414 

15.507 

0.0 

0.0 

-0.235 

0.435 

-0.393 

-0.448 

4.710 

5.252 

0.0 

0.0 

-1.157 

0.743 

0.0 

0.0 

0.840 

0.497 

44.312 

44.236 

0.0 

0.0 

0.662 

1.582 

0.0 

0.0 

0.332 

0.611 

1.202 

1.039 

0.0 

0.0 

15.287 

8.057 

0.0 

0.0 

-2.611 

-0.380 

10.539 

4.325 

0.0 

0.0 

6.708 

9.283 

0.0 

0.0 

73.573 

95.293 

0.0 

0.0 

5.567 

3.686 

0.0 

0.0 

3.611 

2.911 

2.579 

2.771 

7.449 

7.763 

0.0 

0.0 

2.2076 

2.2076 

1.2320 

1.2320 

2.2076 

2.2076 

nergies  are  given  in  ev  ana  aistances  in  Angstroms. 


Figure  9.  Contour  Plot  of  the  N02  2 A'  Potential  Energy  Surface  at  a  NOO  Bond  Angle 
of  110°.  The  Contour  Values  in  Ev  Are  3(1),  2(2),  1(3),  0.3(4),  0.0(5),-0.5(6), 
and  -1(7)  Relative  to  the  N(4S)  +  02  Asymptote. 


Figure  10.  Contour  Plot  of  the  N02  4A'  Potential  Energy  Surface  at  a  NOO  Bond  Angle 
of  110°.  The  Contour  Values  in  eV  are  3(1),  2(2),  1(3),  0.5(4),  0.0(5), -0.5(6), 
and  -1(7)  Relative  to  the  N(4S)  +  02  Asymptote. 
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N-0-0  ANGLE  (deg) 

Figure  11.  Bending  Potentials  of  the  2 A'  and  4 A'  Surfaces  in  the  Region  of  the  Saddle 
Point.  The  Solid  Lines  are  the  Analytical  Fits  Using  Equation  (1)  and  Filled 
Symbols  are  the  Adjusted  ab  initio  Points. 

3.3  Classical  Trajectory  Methodology 

The  quasiclassical  trajectory  (QCT)  method* 54)  is  used  to  compute  the  thermal  reaction  rate 
constants,  and  the  translational  energy  dependence  of  the  reaction  cross  sections  and  the  final 
vibrational/rotational  state  distributions.  The  validity  of  classical  mechanics  for  calculating 
reaction  rate  constants  and  product  distributions  has  been  discussed  elsewhere.*54*  Classical 
calculations  are  expected  to  be  reliable  at  energies  greater  than  the  barrier  heights,  where  quantum 
effects  such  as  tunneling  may  be  important,  for  state-to-state  rate  constants  which  are  strongly 
classically  allowed.  Calculations  are  carried  out  separately  for  the  2 A'  and  4A'  potential  energy 
surfaces  and  then  statistically  averaged  using  the  electronic  degeneracy  factors  of  1/6  and  1/3, 
respectively.  Standard  Monte  Carlo  techniques  are  used  to  compute  the  thermal  reaction  rate 
constant  as  a  function  of  temperature,  where  the  translational  energy,  ET,  and  the  initial 
vibrational-rotational  states,  (v,  j),  are  chosen  from  Boltzmann  distributions.  Cross  sections  and 
final  NO(v',  j')  distributions  are  also  calculated  as  a  function  of  initial  translational  energy  from 
0.3  eV  to  2.75  eV  with  the  initial  (v,  j)  states  selected  from  a  300  K  Boltzmann  distribution.  The 
effect  of  electronic  angular  momentum  on  the  NO  rotational  state  has  been  ignored,  as  the  NO 
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diatomic  properties  appropriate  to  the  NO(X2II1/2)  ground  state  have  been  used.  The  classical 
equations  of  motion  are  integrated  using  a  variable  step  size  predictor-corrector  method.(55)  The 
final  vibrational  and  rotational  quantum  numbers  (v\  j')  are  obtained  from  the  correspondence 
rules  v'd  =  (v'  +  ‘/2)h  and  j'c,  =  (j'  +  where  v'c|  and  j'cI  are  the  classical  vibrational  and 
rotational  action  variables.  The  final  (v',  j')  distributions  are  then  obtained  using  the  standard 
histogram  method.  A  total  of  790,000  (540,000  for  the  2A'  PES  and  250,000  for  the  4A'  PES) 
trajectories  are  used  in  the  current  study. 

3.4  Results  and  Discussion 

3.4.1  Thermal  rate  constants 

Quasiclassical  trajectory  thermal  reaction  rate  constants  are  computed  for  the  temperature 
range  of  500  K  to  5000  K.  A  comparison  of  the  calculated  thermal  rate  constant,  k(T),  with  a 
recommended  fit  to  the  available  kinetic  data(43)  is  shown  in  Figure  12.  The  QCT  thermal  rate 
constant  agrees  well  with  the  experimental  measurements  over  the  entire  temperature  range, 
although  the  calculations  consistently  underpredict  the  rate  constants  below  1250  K.  This 
systematic  disagreement  is  probably  due  to  two  factors:  the  barrier  height  of  0.3  eV  is  either  too 
high,  or  there  is  a  significant  contribution  from  quantum  mechanical  tunneling,  which  cannot  be 
accounted  for  with  the  QCT  method.  Although  important  for  determining  low  temperature 
thermal  rate  constants  and  the  reaction  attributes  at  translational  energies  near  the  barrier,  neither 
one  of  these  concerns  should  be  as  significant  for  calculation  of  the  high  energy  (>0.5  eV)  cross 
sections  and  final  state  distributions  reported  here.  The  temperature  dependence  of  the 
components  of  the  thermal  rate  constant  for  the  2Af  and  4 A'  surfaces  is  shown  in  Figure  13.  As 
expected,  the  contribution  of  2A'  is  dominant  at  temperatures  less  than  1500  K,  due  to  the  lower 
barrier  (0.3  eV)  for  this  surface.  The  4A'  PES  begins  to  make  an  important  contribution  at  2000 
K,  and  is  as  important  as  the  2  A'  surface  at  5000  K.  This  figure  illustrates  the  uncertainty  in 
extrapolating  low  temperature  rate  constants  to  higher  temperature,  where  a  significant 
contribution  to  the  rate  constant  may  come  from  a  mechanism  or  potential  energy  surface  not 
important  at  lower  temperatures.  In  particular,  the  present  calculations  indicate  that  if  the  lower 
temperature  rate  constant  (only  2 A!)  is  extrapolated  to  estimate  the  reaction  rate  constant  at  5000 
K,  the  total  rate  constant  would  be  underestimated  by  a  factor  of  2.  The  present  results  are  also 


28 


in  good  agreement  with  the  rate  constants  calculated  by  Jaffe  et  al.(48^  using  a  completely 
different  analytical  form  for  the  potential  energy  surfaces. 


u 


Figure  12.  Comparison  of  the  Experimental  ( - )  and  QCT(— •— )  Thermal  Reaction  Rate 

Constants  as  a  Function  of  Temperature. 


Figure  13. 


Contributions  of  the  2  A'  and  4 A'  Potential  Energy  Surfaces  to  the  Thermal  Rate 
Constant. 
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3.4.2  Energy  Dependent  Cross  Sections 

A  more  useful  quantity  for  the  thermospheric  NO  production  models* l0,  35'37)  is  the  reaction 
cross  section  as  a  function  of  initial  translational  energy,  ET.  Figure  14  shows  the  total  reaction 
cross  section,  oR(ET),  calculated  as  a  function  of  ET.  Although  the  initial  02  vibrational  and 
rotational  states  were  averaged  over  a  Boltzmann  distribution  at  300  K,  the  cross  sections  are 
insensitive  to  initial  rotational  excitation.  Also  shown  in  Figure  14  is  the  individual  contribution 
of  the  2A'  and  4A'  potential  energy  surfaces  to  the  cross  section.  Note  that  the  statistical 
electronic  degeneracy  factors  have  been  included  in  the  individual  cross  sections.  As  expected, 
the  reaction  cross  section  is  dominated  by  the  2 A'  PES  for  energies  below  1  eV.  As  the  initial 
translational  energy  is  increased,  the  relative  importance  of  the  4A'  reaction  increases  and 
becomes  greater  than  the  2 A!  reaction  above  1 .25  eV  due  to  the  larger  degeneracy  factor. 

3.4.3  Product  Energy  Deposition 

To  assess  the  possibility  of  the  N(4S)  +  02  reaction  contributing  to  the  dayglow  NO  emission 
from  CIRRIS  1A(38,  39)  and  the  NO  emission  from  the  EXCEDE  artificial  auroral  experiment,*41* 
it  is  important  to  establish  the  NO  vibrational-rotational  distributions.  The  calculated  NO 
vibrational  distribution  at  initial  translational  energies  of  0.5  eV  and  1.5  eV  is  presented  in  Figure 
15.  At  the  lowest  translational  energy  (0.5  eV),  the  vibrational  distribution  is  relatively  flat  from 
v'  =  2  to  v'  =  6  and  then  rapidly  decreases  until  the  maximum  vibrational  state  (v'  =  9)  is 
reached.  As  ET  is  increased  to  1.5  eV,  the  maximum  vibrational  state  populated  is  v'  =  14.  The 
vibrational  distribution  peaks  around  v'  =  4  and  then  exponentially  decreases  until  ~v'  =  12.  The 
contribution  of  each  PES  to  the  final  vibrational  distribution  at  1.5  eV  is  shown  in  Figure  16. 
Although  the  vibrational  distributions  for  each  PES  are  normalized  to  unity,  the  reaction  cross 
sections  are  approximately  equal  at  1.5  eV  (Figure  14)  and  the  resultant  distribution  in  Figure  15 
is  nearly  an  average  of  the  two  distributions.  The  greater  extent  of  vibrational  excitation  on  the 
2 A'  surface  than  on  the  4A'  surface  is  consistent  with  the  energy  release*52*  and  curvature  of  the 
minimum  energy  path*53*  of  the  respective  surfaces  as  discussed  earlier.  The  average  vibrational 
quantum  numbers  at  ET=  1.5  eV  are  4.6  and  3.4  on  the  2A'  and  4A'  surfaces,  respectively. 
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Figure  14. 


Figure  15. 
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Contribution  of  the  QCT  Cross  Sections  for  the  2 A'  and  4A'  Potential  Energy 
Surfaces  to  the  Total  Cross  Section  as  a  Function  of  the  Initial  Relative 
Translational  Energy. 
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NO  Vibrational  Distribution  as  a  Function  of  the  Final  Vibrational  Quantum 
Number  for  Initial  Relative  Translational  Energies  of  0.5  and  1.5  eV. 
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Figure  16.  Contribution  of  the  2  A'  and  4 A'  Potential  Energy  Surfaces  to  the  NO 

Vibrational  Distribution  as  a  Function  of  the  Final  Vibrational  Quantum 
Number  at  an  Initial  Relative  Translational  Energy  of  1.5  eV. 


The  present  results  for  the  2 A'  PES  are  in  qualitative  agreement  with  the  calculations  of 
Gilibert  et  al./40^  The  largest  difference  between  the  calculations  is  our  prediction  that  the 
reaction  threshold  occurs  at  ~0.3  eV,  about  0.2  eV  lower  than  the  previous  calculation.  There 
are  also  small  differences  in  the  NO  final  vibrational  distribution,  as  the  present  calculations  show 
that  the  NO  is  formed  with  more  energy  in  vibration.  This  result  is  due  to  the  present  A'  PES 
being  more  attractive  than  the  surface  of  Gilibert  et  al.,  which  is  consistent  with  the  ab  initio 
calculations. 

There  have  been  three  previous  measurements  of  the  NO  vibrational  distribution  at  room 
temperature/44'46^  The  measurements  of  Rahbee  and  Gibson*44^  fall  off  exponentially  with 
increasing  NO  vibrational  quantum  number  in  contrast  with  the  ET  =  0.5  eV  calculation.  Herm 
et  al./45*  have  pointed  out  the  importance  of  considering  the  effects  of  vibrational  quenching. 
The  experiment  of  Herm  et  al./45)  shows  a  NO  vibrational  distribution  that  is  qualitatively  similar 
to  the  QCT  distribution  at  ET  =  0.5  eV,  although  the  experiment  indicates  a  substantial  population 
in  the  v'  =  0  state.  However,  again  due  to  vibrational  relaxation,  it  is  claimed  that  the  measured 
distribution  only  represents  a  lower  bound  to  the  true  nascent  distribution.  An  oscillating  NO 
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vibrational  distribution  was  determined  in  the  most  recent  experiment,*46'  which  disagrees  with 
previous  measurements.  Gilibert  et  al.,*49'  have  attempted  to  explain  qualitatively  the  oscillations 
in  the  vibrational  distribution  by  considering  a  Franck-Condon  model,  although  their  analysis  is 
not  conclusive.  The  present  QCT  calculations  at  0.5  eV  indicate  that  almost  50%  of  the  available 
energy  goes  into  vibration,  somewhat  greater  than  the  experimental  observation  of  24-34%.  Thus, 
the  QCT  calculations  are  consistent  with  the  assertion  of  Herm  et  al.,*45'  that  previous 
experiments  actually  measured  a  partially  relaxed  vibrational  distribution. 

Analysis  of  data  from  the  CIRRIS  1A(38'39'  and  EXCEDE  III(4,)  observations  have  shown  that 
NO  is  formed  in  the  atmosphere  with  extensive  rotational  excitation.  If  the  "hot"  N(4S)  +  02 
reaction  is  important  in  NO  formation  as  indicated  by  Sharma  et  al.,*9'  the  QCT  calculations 
should  show  high  rotational  excitation.  The  NO  rotational  distributions,  summed  over  all  final 
vibrational  states,  are  shown  in  Figure  17  for  translational  energies  of  0.5  eV  and  1.5  eV.  These 
distributions  are  Boltzmann-like  with  rotational  temperatures  of  4000  K  and  10000  K  at 
translational  energies  of  0.5  eV  and  1.5  eV,  respectively.  The  large  amount  of  NO  rotational 
excitation  from  the  N(4S)  +  02  reaction  is  in  quantitative  agreement  with  the  field  experiments. 
The  rotational  distributions  for  v'  of  1,  4,  and  7  are  shown  in  Figure  18,  and  correspond  to 
rotational  temperatures  of  approximately  16000  K,  9400  K,  and  5900  K,  respectively.  The 
average  rotational  quantum  number  at  0.5  eV  and  1 .5  eV  is  shown  as  a  function  of  vibrational 
quantum  number  in  Figure  19.  The  decrease  in  rotational  excitation  as  the  NO  vibrational  state 
increases  is  similar  to  that  expected  from  a  microcanonical  rotational  distribution.  Extensive 
rotational  excitation  of  NO  has  also  been  observed  recently  by  Guadagnini  et  al.,*56)  in  a  detailed 
study  of  NO  formation  from  a  family  of  PINO  reactions. 

Lastly,  the  energy  dependence  of  the  energy  deposition  from  the  N(4S)  +  02  reaction  is 
presented  in  Figure  20  as  the  fraction  of  the  total  energy  available.  At  low  translational  energies, 
the  reaction  deposits  a  large  fraction  (-60%)  of  the  available  energy  in  vibration,  in  agreement 
with  many  studies  of  exothermic  reactions.  As  the  initial  energy  increases  (and  the  initial  orbital 
angular  momentum),  the  incremental  increase  in  translational  energy  is  preferentially  converted 
into  NO  rotational  energy,  which  has  been  observed,  for  example,  by  Muckerman*57'  in  the  F  + 
HD  reaction.  This  behavior  can  be  qualitatively  explained  by  the  spectator  stripping  model. 
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Figure  17.  NO  Rotational  Distribution  as  a  Function  of  the  Final  Rotational  Quantum 
Number  for  Initial  Relative  Translational  Energies  of  0.5  and  1.5  eV.  The 
Solid  Lines  Correspond  to  Boltzmann  Distributions  at  Rotational  Temperatures 
of  4000  K  and  10000  K. 
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Figure  18.  NO  Rotational  Distribution  as  a  Function  of  Final  Rotational  Quantum  Number 
for  Several  Final  Vibrational  Quantum  Numbers  at  an  Initial  Relative 
Translational  Energy  of  1.5  eV. 
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Figure  19.  The  Average  Rotational  Quantum  Number  as  a  Function  of  the  Final 

Vibrational  Quantum  Number  at  Initial  Translational  Energies  of  0.5  and  1.5 
eV. 
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Figure  20.  The  Fraction  of  Total  Energy  in  Translation,  Vibration,  and  Rotation  for  the 
Reaction  Products  as  a  Function  of  Initial  Relative  Translational  Energy. 
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3.5  Summary 


Extensive  quasiclassical  trajectory  calculations  have  been  performed  for  the  N(4S)  +  02(X3Eg) 
— ►  NO(X2II)  +  0(3P)  reaction  using  realistic  ab  initio  potential  energy  surfaces.  The  ab  initio 
calculations  adequately  characterized  the  saddle  point  region  on  both  the  2 A'  and  4A'  electron 
surfaces,  which  is  important  for  the  determination  of  the  thermal  rate  constants  and  energy 
dependent  cross  sections.  However,  the  energy  release  along  the  minimum  energy  path,  which 
strongly  influences  the  final  NO  vibrational  distribution,  was  only  established  for  the  A'  surface. 
The  reaction  attributes  were  obtained  independently  on  the  2A'  and  4A'  electronic  states  and 
statistically  averaged  using  the  respective  degeneracy  factors.  Good  agreement  has  been  obtained 
with  experimental  thermal  rate  constants,  suggesting  that  the  calculations  provide  a  good 
representation  of  this  reaction.  Reaction  cross  sections  and  internal  NO  vibrational-rotational 
distributions  have  also  been  calculated  as  a  function  of  initial  translational  energy.  The  calculated 
NO  vibrational  distributions  at  low  translational  energy  are  consistent  with  available  experimental 
measurements,  although  the  experimental  results  are  colder  than  the  calculated  distribution.  It 
has  been  conjectured  that  the  experimental  distributions  are  partially  relaxed.  At  translational 
energies  greater  than  1  eV,  the  NO  vibrational  distributions  are  uncertain  due  to  the  increasing 
influence  of  the  4A'  surface  on  the  energy  disposal.  Recent  field  measurement  programs 
observed  NO  distributions  with  extensive  vibrational  excitation  characterized  by  rotational 
temperatures  on  the  order  of  4000  K  to  7000  K,(39,41)  which  the  present  QCT  calculations  show 
are  consistent  with  NO  formation  from  the  N(4S)  +  02  reaction.  Further  calculations  of  NO 
formation  in  the  atmosphere  from  translationally  hot  N(4S)  +  02  using  the  present  results  should 
provide  more  quantitative  evaluation  of  the  relative  importance  of  N(4S)  as  a  source  of  NO  in 
the  thermosphere. 
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4.  EVIDENCE  FOR  NO  FORMATION  FROM  HYPERTHERMAL  N(4S) 
ATOMS  IN  AN  ARTIFICIAL  AURORAL  EXPERIMENT 

4.1  Introduction 

This  section  discusses  the  analysis  of  electron-beam  enhanced  NO(Av  =  1)  emission  observed 
by  the  interferometer  during  EXCEDE  III.  In  addition  to  considering  the  traditional  N(2D)  +  02 
reaction  at  thermal  energies,  the  chemistry  of  translationally  hot  N  atoms  in  NO  formation  has 
been  investigated  in  detail,  with  emphasis  on  the  hyperthermal  N(4S)  +  02  reaction.  It  is  shown 
that  a  quantitative  description  of  rotationally  thermal  and  hot  NO  emission  observed  in  EXCEDE 
III  requires  the  inclusion  of  hyperthermal  N  atom  chemistry  in  our  chemical  kinetics  model. 

Section  2  discusses  the  evidence  for  the  importance  of  the  reactions  of  hyperthermal  N(4S) 
and  N(2D)  with  02  in  producing  NO.  The  chemical  dynamics  of  the  N(4S)  +  02  reaction  is 
discussed  in  Section  3,  and  preliminary  calculations  of  NO  vibrational  populations  using  a 
hyperthermal  nitrogen  atom  kinetics  model  is  presented  in  Section  4.2.  A  brief  summary  is  given 
in  Section  4.3. 

4.2  Chemical  Kinetics  Model 

The  reactions  of  metastable  nitrogen  atoms  with  molecular  oxygen  are  thought  to  be  the  major 
source  of  nitric  oxide  chemiluminescence  in  the  thermosphere/51  Earlier  models  of  NO  formation 
assumed  that  the  chemistry  of  N(2D)  is  the  most  significant  contributor  to  NO  formation,  while 
the  relaxation  of  N(2P)  to  N(2D)  and  the  reaction  of  N(2P)  with  02161  make  a  relatively  minor 
contribution  in  the  thermospheric  dayglow1331  and  aurora/14-  32) 

Reaction  (1)  has  previously  been  considered  as  a  source  of  rotationally  thermal  NO/7,  81  The 
production  mechanism  of  rotationally  hot  NO  is  currently  unknown,  although  it  has  been 
suggested  that  the  reaction  of  N(2P)  with  02  is  a  possible  source/81  However,  it  has  been  shown 
(J.  W.  Duff,  unpublished  work,  1994)  that  Reaction  (2)  can  not  be  a  significant  source  of 
rotationally  thermal  or  hot  NO  emission  under  the  conditions  of  the  EXCEDE  III  experiment. 

Further  evidence  of  the  inadequacy  of  Reaction  (1)  is  provided  by  the  analysis  of  the 
interferometer  data  at  apogee  which  shows  significant  emission  from  highly  rotationally  excited 


37 


NO.  At  apogee,  Reaction  (1)  does  not  produce  the  correct  NO  vibrational  distribution  to  explain 
the  rotationally  thermal  populations  obtained  from  the  interferometer  data.  Recently,  it  has  been 
suggested  that  the  reaction  of  translationally  hot  N*(4S)  with  02  may  be  an  important  contributor 
to  NO  formation  in  the  thermosphere/34, 36)  Furthermore/9^  have  argued  the  importance  of  the 
N*(4S)  +  02  reaction  in  recent  observations  of  highly  rotationally  excited  NO  vibrational 
emissions  in  the  dayglow/38,  39^ 

In  order  to  access  the  importance  of  Reaction  (3)  as  a  significant  source  of  NO  which  is  both 
highly  rotationally  and  vibrationally  excited  in  EXCEDE  III,  we  have  created  a  time-dependent 
kinetics  model  incorporating  a  mechanism  for  producing  NO  from  "hot"  N  (4S).  The  most 
important  inputs  for  this  model  are  the  reaction  cross  section  for  the  N*(4S)  +  02  — ♦  NO(v,j)  + 
O  process,  and  the  quenching  of  N*(4S)  by  major  atmospheric  species.  The  model  used  for  the 
quenching  of  N*(4S)  is  similar  to  that  used  by  Logan  and  McElroy  in  the  study  of  fast  oxygen 
atoms/58*  The  translational  energy  distribution  function  for  the  N  (4S)  +  02  reaction  is  defined 
in  the  energy  range  0-1.5  eV.  This  energy  range  is  subdivided  into  finite  bins  of  width  0.1  eV. 
The  reaction  rate  constants  for  N*(4S)  +  02  from  the  quasiclassical  trajectory  calculations  are  then 
computed  at  a  translational  energy  corresponding  to  the  average  energy  of  each  bin.  The  various 
sources,  and  associated  translational  energy,  of  N*(4S)  from  irradiation  of  the  atmosphere  are 


e  +  N  -  e  +  N(4S)  +  N(2D,2P)  [0  -  2  eV]  (28) 

pip 

n!  +  e  ~  aN(*S)  +  (2  -  a)N(2D)  [1.7  eV\  (29) 

O*  +  N2  -  NO*  +  N(4S)  [0.75  eV\  (30) 

N+  +  0  -  O  +  N(4S)  [1.7  eVJ  (31) 
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[1.5  eVJ 


(32) 


NO+  +  e  -  O  +  N(4S ) 

S 


O  +  N(2D)  -  O  +  W(4S)  [1.27  eV\  (33) 


N*  +  O  -  O*  +  N(*S)  [0.5  eVJ  (34) 

O  +  N(2/>)  -  O  +  A^(45)  [1.9  eF]  (35) 

The  translational  energies  from  Reactions  (28)  -  (35)  are  transformed  to  the  N(4S)  +  02  center 
of  mass  and  placed  in  the  appropriate  energy  bin.  There  are  estimates  of  the  cross  section  for 
the  quenching  of  N*(4S),  although  the  uncertainties  are  unknown.  We  have  chosen  a  reasonable 
value  of  5  x  10‘16  cm2  for 


N(4S)  +  M  -  N(4S)  +  M  (36) 

which  is  to  be  compared  with  the  gas  kinetic  estimate  of  3  x  10'15  cm2. 

The  rotationally  hot  vibrational  populations  calculated  from  Reaction  (3)  are  compared  to  the 
NO  populations  obtained  from  the  interferometer  data*59)  at  an  altitude  of  103  km  (upleg)  in 
Figure  21.  The  excellent  agreement  between  the  chemical  kinetics  model  and  data  provide  the 
first  quantitative  evidence  that  Reaction  (3)  is  the  source  of  rotationally  hot  NO  in  disturbed 
atmospheres.  Also  shown  in  Figure  21  is  a  comparison  of  the  vibrational  populations  obtained 
from  Reaction  (1)  using  the  nascent  vibrational  distribution  of  Kennealy  et  al./7^  with  the 
rotationally  thermal  populations  from  the  interferometer.  Again  the  agreement  of  the  kinetics 
model  and  the  data  is  excellent,  indicating  that  the  N(  D)  is  thermalized  at  103  km  (upleg). 
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A  similar  comparison  of  the  rotationally  thermal  and  hot  NO  vibrational  populations  from  the 
model  and  interferometer  data  at  115  km  (apogee)  is  shown  in  Figure  22.  For  the  thermal 
rotational  component,  the  vibrational  populations  from  the  kinetic  model  are  approximately  a 
factor  of  3  below  the  data.  This  comparison  implies  that  the  N(2D)  is  probably  not  thermalized 
at  115  km,  and  the  rate  constant  and  vibrational  distribution  used  for  Reaction  (1)  is  not 
correct.  Figure  22  also  contains  a  comparison  of  the  rotationally  hot  vibrational  populations  from 
Reaction  (3)  with  the  interferometer  data.  Although  the  kinetic  model  is  in  reasonable  agreement 
with  the  data,  the  model  underpredicts  the  vibrational  populations  for  low  vibrational  quantum 
numbers.  At  this  altitude  (115  km),  N(4S)  translational  energies  in  the  range  of  1  to  2  eV  is  most 
important  for  producing  NO.  As  is  evident  from  Figure  14  in  Section  3,  the  4A'  PES  makes  a 
significant  contribution  to  the  NO  vibrational  distribution.  Unfortunately,  the  regions  of  the  PES 
which  are  most  important  for  determining  the  vibrational  distribution  on  the  4A'  PES  are  not 
available  from  the  ab  initio  calculations.  Evidently  the  4A'  PES  is  more  repulsive  (i.e.,  colder  NO 
vibrational  distribution)  than  the  PES  used  in  Section  3.  Testing  of  this  hypothesis  will  have  to 
await  more  extensive  ab  initio  calculations. 
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Figure  21.  Comparison  of  the  NO  Vibrational  Populations  for  the  Rotationally  hot 

Component  Derived  from  the  Interferometer  Data  at  an  Altitude  of  103  km 
During  Upleg  with  the  EXCEDE  Kinetic  Model  Using  the  N*(4S)  +  02 
Mechanism. 


0  123456789  10  11  12 


VIBRATIONAL  QUANTUM  NUMBER 

Figure  22.  Same  as  Figure  21  Except  for  an  Altitude  of  115  km  (apogee). 
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4.3  Summary 


This  analysis  has  provided  the  first  quantitative  evidence  of  the  importance  of  hyperthermal 
N(4S)  and  N(2D)  atoms  in  the  formation  of  vibrationally  and  rotationally  excited  NO.  The  rate 
constants  and  vibrational  distributions  for  rotationally  hot  NO  are  obtained  from  extensive 
quasiclassical  trajectory  calculations  for  the  N(4S)  +  02(X3Eg)  — ♦  NO(X2II)  +  0(3P)  reaction 
using  realistic  ab  initio  potential  energy  surfaces.  Good  agreement  has  been  obtained  with 
experimental  thermal  rate  constants,  indicating  that  the  calculations  are  a  good  representation  of 
this  reaction.  Reaction  rate  constants  and  internal  NO  vibrational -rotational  distributions  have 
also  been  calculated  as  a  function  of  initial  translational  energy.  The  NO  vibrational-rotational 
distributions  at  low  translational  energy  are  consistent  with  available  experimental  measurements. 
Recent  field  measurement  programs  observed  NO  distributions  with  extensive  vibrational 
excitation  characterized  by  rotational  temperatures  on  the  order  of  4000  K  to  7000  K./39'  41  ^  which 
the  present  QCT  calculations  show  are  consistent  with  NO  formation  from  the  N(4S)  +  02 
reaction. 

Quantitative  agreement  between  the  chemical  kinetics  model  developed  for  EXCEDE  and  the 
vibrational  populations  derived  from  the  interferometer  data  is  obtained  under  conditions  of 
thermalization  of  nitrogen  atoms  (i.e.,  at  103  km  under  max  dose  conditions).  Analysis  of  the 
vibrational  populations  from  the  interferometer  under  other  conditions  indicate  that  either 
hyperthermal  N(2D)  atoms  also  play  an  important  role  in  NO  formation  and/or  additional 
mechanism(s)  are  important.  Information  concerning  the  N(2D)  +  02  reaction  at  hyperthermal 
translational  energies,  analogous  to  that  obtained  in  the  present  work  for  the  N(4S)  +  02  reaction, 
is  not  currently  available. 

The  results  of  this  study  indicate  that  a  quantitative  understanding  of  NO  formation  in  the 
thermosphere  will  require  a  detailed  investigation  into  the  dynamics  of  N(  D)  +  02  reaction. 
Limited  ab  initio  studies  of  the  N(2D)  +  02  system  have  been  carried  out  to  characterize  the 
reaction  pathways.(60)  Although  a  total  of  five  potential  energy  surfaces  will  have  to  be 
considered  to  treat  the  dynamics  of  N(2D)  +  02,  the  calculations  are  well  within  the  scope  of  the 
quasiclassical  trajectory  method  discussed  in  this  report.  Given  a  more  extensive  and  well  chosen 
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set  of  ab  initio  points  for  the  PES  than  is  currently  available,  such  calculations  for  this  system 
could  be  pursued  in  the  future. 
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5.  SUMMARY 


The  EXCEDE  III  experiment  has  provided  quantitative  characterization  of  the  N(2D)  5200 
A,  the  N(2P)  3466  A,  and  the  NO(Av=l  and  Av=2)  emissions  resulting  from  the  irradiation  of 
the  lower  thermosphere  (90  -  1 15  km)  by  an  intense  electron  beam.  A  detailed  chemical  kinetics 
model  has  been  developed  describing  the  interaction  of  an  electron  beam  with  the  atmosphere 
including  the  geometrry  of  the  EXCEDE  III  measurement.  The  model  indicates  that  the  N(  D) 
emission  can  be  quantitatively  explained  by  assuming  accepted  production  rates  of  N(2D)  and  N2 
created  by  the  electron  beam  dissociation,  consistent  with  previous  auroral  models.  Using  the 
previous  estimate  of  0.22  N(2P)  per  ion  pair(,4)  from  electron  dissociation  of  N2  results  in  a 
model  prediction  which  is  approximately  50%  below  the  data,  in  disagreement  with  previous 
analyses  of  N(2P)  auroral  and  dayglow  emission.* I4,32,33^  However,  it  should  be  recalled  that  the 
EXCEDE  III  measurement  of  N(2P)  was  not  made  under  steady  state  conditions,  which  are 
appropriate  to  the  atmospheric  observations.  The  present  analysis  would  imply  that  the  amount 
of  N(2P)  from  N2  dissociation  should  be  about  a  factor  of  2  greater  than  previously  assumed  or 
an  additional  source  of  N(2P)  exists  for  intense  electron  beams  such  as  that  used  in  EXCEDE  III. 

To  assess  the  contribution  of  the  nonthermal  N(4S)  to  NO  production,  extensive  quasiclassical 
trajectory  calculations  have  been  performed  for  the  N(4S)  +  02(X3S'g)  — *•  NO(X2II)  +  0(3P) 
reaction  using  realistic  ab  initio  potential  energy  surfaces.  Good  agreement  has  been  obtained 
with  experimental  thermal  rate  constants,  suggesting  that  the  calculations  provide  a  good 
representation  of  this  reaction.  Reaction  cross  sections  and  internal  NO  vibrational-rotational 
distributions  have  also  been  calculated  as  a  function  of  initial  translational  energy.  Recent  field 
measurement  programs  observed  NO  distributions  with  extensive  vibrational  excitation 
characterized  by  rotational  temperatures  on  the  order  of  4000  K  to  7000  K,*39'41*  which  the 
present  QCT  calculations  show  are  consistent  with  NO  formation  from  the  N(4S)  +  02  reaction. 

The  first  evidence  of  the  importance  of  hyperthermal  N(4S)  and  N(2D)  atoms  in  the  formation 
of  vibrationally  and  rotationally  excited  NO  is  described  in  this  report.  The  rate  constants  and 
vibrational  distributions  for  rotationally  hot  NO  are  obtained  from  the  quasiclassical  trajectory 
calculations  for  the  N(4S)  +  02(X3Sg)  discussed  in  Section  3.  Quantitative  agreement  between 
the  chemical  kinetics  model  and  the  vibrational  populations  derived  from  the  interferometer  data 
is  obtained  under  conditions  of  thermalization  of  nitrogen  atoms  (i.e.,  at  103  km  under  max  dose 
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conditions).  Analysis  of  the  vibrational  populations  from  the  interferometer  under  other 
conditions  at  apogee  indicate  that  either  hyperthermal  N(2D)  atoms  also  play  an  important  role 
in  NO  formation  and/or  additional  mechanism(s)  need  to  be  considered. 
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